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Abstract 

The influence of relativity on the triton binding energy is investi- 
gated. The relativistic three-dimensional version of the Bethe-Salpeter 
equation proposed by Blankenbecler and Sugar (BbS) is used. Rel- 
ativistic (non-separable) one-boson-exchange potentials (constructed 
O 

rH \ in the BbS framework) are employed for the two-nucleon interaction. 

•'^ . In a 34-channel Faddeev calculation, it is found that relativistic ef- 

X 

C^ ' fects increase the triton binding energy by about 0.2 MeV. Including 

charge-dependence (besides relativity), the final triton binding energy 
predictions are 8.33 and 8.16 MeV for the Bonn A and B potential, 
respectively. 
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1 Introduction 

One of the main objectives of nuclear few-body physics is to learn something 
about nuclear forces. When the few-body problem can be solved accurately, 
results can be related directly to the elementary (two- and many-body) forces 
used as input in the calculations. Since the two-nucleon scattering data 
determine the nucleon-nucleon (NN) interaction only on-shell, one has to 
consider A-nucleon systems with A > 2 to test off-shell properties of nuclear 
potentials. Similarly, many-body forces can only be tested in systems with 
A > 2. Guided by these ideas, comprehensive work has been performed 
for, particularly, the trinucleon (3A^) system. For this system, an exact and 
practically feasible theory has been developed by Faddeev ^. 

In recent years, many benchmark calculations of the triton binding energy 
employing realistic A^A^ potentials have been performed p|, ^, Q — with the 
objective of obtaining (indirect) information on the off-shell behaviour of the 
two-nucleon potential and to assess the size of the nuclear three-body force. 
It has been shown that when two-nucleon forces only are applied, the triton 
is underbound by an amount which varies from 0.2 to 1.1 MeV (depending 
on which NN potential is used in the calculation). The conventional interpre- 



tation of this result is that it 'proves empirically' the existence of three-body 
forces (which would contribute 0.2 MeV or more to the triton binding). As 
shown in many calculations f^ , it is then possible to patch up the discrepancy 
with experiment by including phenomenological three-body forces which can 
easily be fit to any amount of missing binding energy. Unfortunately, not 
much physics can be learned in this ad hoc fashion. 

Strictly speaking, most many-body forces are an artefact of theory. They 
are created by freezing out degrees of freedom contained in the full Hamil- 
tonian of the problem under consideration. This fact suggests how one may 
learn some physics from three-body-force calculations: when a new degree 
of freedom is introduced (besides the nucleon), it should be taken into ac- 
count consistently in the two- and the three-body problem. The theory of 
the NN interaction has shown that isobar degrees of freedom (particularly 
the A(1232)) are crucial for a realistic and reliable description of the nuclear 
force . Trinucleon calculations taking the A-degree of freedom consistently 
into account have been performed first by the Hannover group [0. Recently, 
this type of calculation has been updated and improved pf. The new re- 
sult is that the sum of all three-body force contributions arising from the 
A-isobar is repulsive; the repulsive contribution to the triton binding energy 



varies between 0.06 and 0.4 MeV depending on the A model [Q. This small 
total result is due to a cancelation between an attractive three-body force 
contribution of about 1.5 MeV and a slightly larger repulsive dispersive ef- 
fect on the A-diagrams contained in the two-nucleon interaction. Even if the 
7T — N S-wave part of a 27r-exchange three-body force |^ (not included in a A 
model) is added, the total result is essentially vanishing. Thus, there is little 
hope that a consistent theory for two- and three-body forces will explain the 
triton binding energy within the conventional framework. 

On the background of the above discussion, it becomes increasingly im- 
portant to investigate aspects left out in conventional three-nucleon calcu- 
lations. One of these aspects is relativity. Most present calculations using 
realistic two-body forces are performed in a non-relativistic framework. Thus, 
one may ask whether relativistic effects are responsible for the missing bind- 
ing. Since the triton wave functions contain high momentum components, 
it is by no means clear that the non-relativistic approximation is justified. 
Furthermore, some non-relativistic predictions for the triton binding energy 
(obtained with two-body forces only) are only 0.2 MeV off the empirical 
value [0. Thus, even if relativistic effects are small, they may be large as 
compared to the small gap to be filled by 'non-conventional' nuclear physics. 



In summary, the current status of the triton binding energy program calls 
for an accurate knowledge of the relativistic corrections. 

In the past 25 years, there have been numerous efforts to estimate rel- 
ativistic effects in the three-nucleon system. Roughly, one can distinguish 
between two lines of research. In a model independent approach, one starts 
from the Poincare algebra for the ten generators of the Poincare group, ex- 
pressed in terms of the degrees of freedom of n particles only. This frame 
work is applied in Refs. [10-13] and critically examined in Refs. [14-16]. 

An alternative approach is to use a field-theoretic framework and the 



Bethe-Salpeter equation |T7] or one of its three-dimensional reductions. Pi- 



oneering work along this line has been performed by Tjon and cowork- 



ers 



18| , p!9| . In reference ||I9[ , a Blankenbecler-Sugar ||2^ three-body equation 
for a modified Reid potential is solved. In a series of recent studies, Rupp 
and Tjon solved the full four-dimensional three-body Bethe-Salpeter equa- 



tion |2l|, |22[. However, to make this tremendous numerical project feasible, 
finite-rank separable NN potentials had to be used. In these studies, rela- 
tivistic effects are found to increase the triton binding by about 0.3 MeV. 

In the present paper, we investigate the influence of relativity in the 
3A^ bound- state system using the relativistic meson-theoretic Bonn poten- 



tial p, ^. It is important that relativistic few-body calculations are per- 
formed consistently. By this we mean that the same relativistic formalism has 
to be used for the two- and the three-body systems. Since the relativistic po- 
tential we apply is constructed in the framework of the Blankenbecler-Sugar 
equation ||20| (a relativistic three-dimensional reduction of the Bethe-Salpeter 



equation [T^), we use an equivalent relativistic equation for the three-nucleon 
system. 

Various three-dimensional relativistic three-body equations have been 
proposed, which all have in common the elimination of the relative time com- 
ponent in the original four- dimensional integration over internal momenta. 
This is usually done by replacing, in the three-particle Bethe-Salpeter equa- 
tion, the original Green's function by one that involves a delta function. This 
procedure is carried out in such a way as to preserve relativistic invariance, 
three-particle unitarity and, of course, the correct non-relativistic limit. 

As mentioned above, our study is performed within such relativistic three- 
dimensional approaches, combining, however, a non-separable relativistic 
two-body potential (the Bonn potential), with relativistic Faddeev equations. 
Among our objectives it is to understand to which extent, if at all, the type of 
two-body potential used is related to the size of the relativistic effect. Finally, 



our main motivation is that the influence of relativity on the 3A^ bound-state 
properties is a basic and fundamental problem, and we must achieve a solid 
understanding of it. 

In the next Section, we will briefly recall the relevant equations and de- 
scribe the relativistic corrections. In Section 3 we present and discuss our 
results. 

2 Formalism 

The relativistic two-body problem is described by the Bethe-Salpeter (BS) 
[ P^ equation, which can be extended to the three-body case. Similarly to 
the non-relativistic formalism, the relativistic three-body BS equation can 



be cast into Faddeev form P0| , |2J, ^, ^. Unlike in the non-relativistic case. 



however, one obtains a four-dimensional integral equation, (the additional 
integration variable is due to the intermediate particles being off their mass 
shell), which drastically complicates the numerical solution. To avoid this 
difficulty, various three-dimensional reductions have been worked out as a 
practical generalization of the relativistic Faddeev equations. Here, we just 
briefly recall the procedure. The relativistic Faddeev equations are written 



as 

T = U + Y.UGo^T^ (1) 

with i,j = 1, 2, 3. ti are two-body scattering operators describing all interac- 
tions between particles j and k, and Goi is the free two-particle Green's func- 
tion of the pair j, k acting in the three-body space. (In the non-relativistic 
case, it can be shown that Gqi = G^02 = Gq^, all particles being on their mass 
shell.) The relativistic propagator, e.g. for particles j and k, is then 

{k'j —m? + ie) (kl — iv? + it) 

with ki the four-momentum of particle i in the three-body center-of-mass 
(CM.) frame. The usual procedure is then to replace the relativistic four- 
dimensional propagators by a three-dimensional one which is derived from 



a dispersion integral [|T^, |I9[. This amounts to introducing a Green's func- 
tion which produces the proper right-hand (unitarity) cut. The procedure 
is, however, not unique, and many different versions can be found in the 
literature. 

Instead of the Green's function Eq. (2), one may write a dispersion inte- 
gral in terms of the two-particle invariant energy, ctj, (as we choose to do in 



this work), namely 

/•OO ^ly' 

g, = 4 -r^^mj'^ + k,,f - m'] x 5\{P'^J2 - k,,f - m'\ (3) 

J Am? O^ — Gi ■' 

where fej^ = ikj — kk)/2^ o"j = (P — kiY = Pfk with P = ki + kj + k^ the total 
four-momentum, Pjk = kj + kk, and Pj;. = Ja'JaiPjk- In the three-body 
CM. frame, P = (y^, 0) with ^/s = 3m — i^^^ the invariant mass of the 
three-body system {m denotes the mass of a free nucleon and E^ the binding 
energy of the three- nucleon system). Clearly, this choice cannot be unique, 
since any function which does not contribute to the right-hand cut of a can 
be added without affecting unitarity. Evaluation of the integral in Eq. (3) 



results in the three-dimensional propagator [|T8], |^ 



^ ± 1 . . 

^ o;p4a;2 + |q|2-(yi-^J2 U 

Here, q(q') are the momenta of the spectator particle in the initial (final) state. 



respectively, cUp = Jm? + |pp, p = — |q — q', and ujq = Jm? + |qp. Notice 
that this is just a more convenient notation which allows to drop particles 
indices, namely q = kj, q' = k'^, where i and k are understood to be the 
spectator particles in the initial and final state respectively. Clearly, q and 
q' are also the relative momenta between the spectator and the interacting 
pair in the three-body CM., see Fig. 1. 



Alternatively, one could have written the dispersion integral in terms of 
the three-particle invariant energy, that is 

9i = ^l^^ y—J{[iP'-h)/2 + k,,]'-m'}x5{[iP'-h)/2-k,,]'-m'} (5) 



with kjk as above, s = P^, P' = PJs'/s. The resulting three-dimensional 
propagator would then be |TB|, pT 



2 UJ1+U2+UJ3 , , 

3 = 7 , \ ^2 (6) 

UJ1UJ2 [UJl +UJ2+ 003)^ — S 



with tui = Jm? + |q'P, tUg = Jrn? + |q+ q'P and 0^3 = Jm? + |qp. 

Other alternatives have also been suggested. One could write as well the 
dispersion integral in terms of the one-particle invariant energy, or introduce 
relativistic kinematics in the free resolvent of the non-relativistic Faddeev 
equation; in the latter case, phase space factors are introduced to make the 



volume element of the integration relativistically invariant |18|. All choices 
are equivalent in the sense that they satisfy three-body unitarity and preserve 
the correct non-relativistic limit. This non-relativistic Faddeev propagator 
reads in our notation 



9NR - — I„I2 , 3|„I2 , ™Z7 C^) 

4I 



m|pP + ||qP + mEb 



Besides the modifications in the propagator, additional corrections must 



be included to properly define the dynamical variables. Clearly, in the three- 
body calculation, these are defined in the overall CM. system and must 
appear correctly in the input two-body t-matrix. The diagram in Fig. 1 
is a typical graphical representation of the three-particle reaction i{jk) —>■ 
k{ij). The double lines represent interacting pairs (with momenta — q, — q', 
respectively). The total invariant mass of the two-body sub-system is 

a={^s~E,Y-\q\' (8) 



with Eg = Jm? + |qp, so that {^/s — Eg) is nothing but the energy of the 
pair with total three-momentum — q (righthandside of diagram in Fig. 1). 

Also, the relative momenta of the interacting pair must be defined in 
a relativistic way (in a non-relativistic theory, they would be the Galilean 
invariant relative momenta). It has been shown shown ||2^, ^ that the 



relativistic relative momentum, (which reduces to the momentum of one of 
the particles in their CM. system), is given by the combination 

Q = q + p(|qUq'U)q' (9) 

(for the lower vertex in Fig. 1). In the last expression, z is the angle between 
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q and q', and the function p is given by [^, ^ 

p(|q|, |q'U) = C/\E, + ^-^^ ) (10) 

with ^g/ = {Eq + -Eq+q') — q' • (A completely analogous definition applies 
to the upper vertex, but with q, q' interchanged.) Notice that Q, Q' are 
formally identical to the Galilean invariant relative momenta, if the function 
p is replaced by a factor 1/2. 

Finally, we mention that the relativistic two-body t-matrix, treh is related 
to the non-relativistic one, t^R, by 



trel = tNRjE/mJE'/m (11) 



with E = \Jvn?' + Q^ and E' = Jm? + Q'^. If in the relativistic Blankenbecler- 
Sugar (BbS) equation, t^ei is replaced by Eq. (11) and the quasi-potential, 
Vrel^ is replaced by 



Vrei = VNRjE/mJE'/m (12) 



then the usual non-relativistic Lippmann-Schwinger equation is obtained [^ . 
In summary, the relativistic calculation of the triton binding energy differs 
from the non-relativistic one in three points: 

1. The relativistic Faddeev propagator, Eq. (4), is used instead of the 
conventional non-relativistic propagator, Eq. (7). 
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2. The kinematical variables of the two interacting nucleons are defined 
in a covariant way (cf. Eqs. (8) and (9)). 

3. An invariant two-body t-matrix is used, cf. Eq. (11). 

It will turn out that the last point (invariance of the t-matrix) has the largest 
quantitative effect and is responsible for the increase of the triton binding 
energy. 

3 Results and Discussion 

In Table 1 we present our results for the non-relativistic and relativistic calcu- 
lations of the triton binding energy taking 5 and 34 three-body channels p^ 



into account. We apply the Bonn A and the Bonn B potential |^ for the 



two-body interaction. The relativistic calculations are based upon the for- 
malism outlined in the previous section, with the three-nucleon propagator 
given in Eq. (4). For the solution of the non-relativistic as well as relativistic 
Faddeev equation, we use standard momentum-space techniques [pT |. 



We find a relativistic correction to the triton binding energy of 0.19 MeV 
for the Bonn A and 0.20 MeV for the Bonn B potential. As indicated above. 
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the largest effect comes from tlie use of an invariant t-matrix; tliis increases 
tlie binding energy by 0.26 MeV. Tlie effect of using covariant kinematical 
variables is very small and further increases the binding energy by 0.07 MeV, 
while the relativistic Faddeev propagator, Eq. (4), reduces the binding energy 
by 0.14 MeV. (The numbers given here for the partial effects refer to the 34- 
channel calculation with Bonn A.) 

Using the Bonn A potential, the non-relativistic 34-channel result of 8.32 
MeV is increased to 8.51 by relativity. These numbers refer to the charge- 
independent calculation in which the neutron-proton potential is used for 
all states. If the well-established charge- dependence in the ^Sq two-nucleon 
potential is taken into account, the triton binding energy is reduced by 0.18 
MeV, yielding the final result of 8.33 MeV for the Bonn A potential. For Bonn 
B, the final result, including relativistic effects and charge-dependence in the 
^So state, is 8.16 MeV (cf. Table 1). In our calculations, charge dependence 



is taken into account as described in Ref. [32]. 

The relativistic effect we find is small but not negligible. Another sub- 
tle effect that is usually neglected in triton binding energy calculations is 
charge- dependence. Accidentally this effect is of the same magnitude as the 
relativistic effect. For this reason, we also consider charge- dependence in 

13 



our calculations. Thus, we can provide rather 'complete' final results for the 
triton binding energy predictions from the two-body force only. This final 
result is that the triton is underbound by 0.15-0.30 MeV depending on the 
choice of the relativistic two-nucleon potential. This means that these po- 
tentials almost completely explain the triton binding energy in terms of just 
the two-nucleon force. 

At this point, one may raise the question whether it is reasonable to ob- 
tain almost all the triton binding energy from the two-body force alone. The 
crucial issue here is what to expect from three-body forces. As mentioned 
in the Introduction, the recent investigations by Picklesimer et al. Q have 
confirmed the earlier findings by the Hannover group 0] that a consistent 
treatment of two- and three-body forces results in an almost vanishing to- 
tal contribution to the triton binding from three-body forces. In this light, 
the triton binding can only be understood as resulting essentially from the 
two-nucleon interaction (unless one wants to invoke subhadronic degrees of 
freedom) . 

Thus, crucial for the understanding of our result is the nature of the 
two-body force. It is well-known |^ that a weaker tensor force component 



in the two-nucleon force implies more attraction in few- and many-body 
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systems. This is why Bonn A predicts more binding than Bonn B; and, 
even more so, this is the reason why both Bonn A and B give substantially 
more binding energy than other conventional potentials which, in general, 
predict only about 7.5 MeV. On the other hand, the tensor force should only 
be as weak as allowed by the requirement that the potential be 'realistic'. 
The usual interpretation of the term 'realistic' is that the potential should 
reproduce the two-nucleon scattering data up to pion-production threshold. 
This is, in fact, true for both Bonn A and B, with the exception of the ei 
mixing parameter at energies larger than 100 MeV, which is insufficiently 



reproduced by Bonn A, but correctly predicted by Bonn B ||33[. On those 
grounds, one may be tempted to cast some doubt on the Bonn A three-body 
result. 

To properly discuss this issue, let us first focus on the Bonn B predic- 
tions. The Bonn B potential is with no doubt realistic (in its predictions of 
all two-body bound state and scattering observables, including the ei), and 
nevertheless predicts substantially more triton binding than other realistic 
potentials. For example, the Paris and the Bonn B potential yield almost 
identical predictions for phase shifts, e parameters, and observables Bl], i. 



e. they are as identical on-shell as two different potentials can be. However, 
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the triton binding energy prediction by Bonn B is about 0.4 MeV larger 
than the one by Paris (the charge-dependent Paris resuh is 7.59 MeV |0). 
This clearly reveals off-shell differences between the two interactions. The 
Paris potential is, apart from a p^ operator, a local potential, while Bonn A 
and B are defined in terms of the relativistic momentum-space Feynman am- 
plitudes for one-boson-exchanges (which are clearly non-local expressions). 
These off-shell differences can already be seen in the D-state probability of 
the deuteron which is 5.8% for Paris and 5.0% for Bonn B. 

Concerning now the Bonn A potential, one has to understand the rel- 
evance of the ei parameter for nuclear ground state predictions. In Refer- 
ence [^ , a family of potentials has been considered that has the same low 



D-state probability of the deuteron as Bonn A, but larger ei parameter at 
intermediate energies. It is found that the triton binding energy is always the 
same and, thus, is not affected by the behaviour of the ei above 100 MeV. 

Finally, some comments are in place concerning the relativistic framework 
we use. Given the fact that our relativistic one-boson-exchange two-nucleon 
potentials are constructed within the BbS scheme, the relativistic three-body 
equation derived within the BbS framework certainly is the most consistent 
choice. 
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On the other hand, for the relativistic three-nucleon problem, (unhke the 
non-relativistic case), predictions from state-of-the-art calculations do not 
corroborate and are, in fact, quite method-dependent. For instance, Rupp 



and Tjon have shown ||2T|, ^ that the BbS approximation yields about 75% 
of the relativistic effect as obtained from the four- dimensional BS equation. 
The reason for this difference may be the effective three-body forces due 
to retardation which are automatically included in the four-dimensional BS 



equation and do not exist in static approaches. I^T], ^ 

Naively, one would believe that the BS equation is a priorithe 'best' equa- 
tion. Based upon this believe, the quality of a relativistic three-dimensional 
equation would then be determined by how close it reproduces results ob- 
tained from the BS equation. 



However, as pointed out by Gross ||3^ , there are physical reasons why this 



seemingly most plausible approach may not be correct. For example, the BS 
equation in ladder approximation (that is, with a kernel which is restricted 
to one-boson-exchanges) does not generate the desired one-body equation in 
the one-body limit (i. e. when one of the particles becomes very massive), 
while a large family of three-dimensional relativistic reductions of the BS 
equation does have the correct one-body limit. If the sum of all crossed 
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ladders is added to the ladders, then the correct limit is obtained in the BS 
equation. As a matter of fact, in almost all applications in nuclear physics 
(including the present work) the kernel is restricted to one-boson-exchanges. 
Furthermore, when going beyond the second order kernel, there is no rapidly 



converging series of irreducible kernels in the BS framework |]36|. Better 
convergence is obtained for some three-dimensional equations. 

Thus, the BS equation may not be the optimal choice for relativistic 
investigations in nuclear physics. Given the problems intrinsic to relativis- 
tic approaches, three-dimensional relativistic equations and results obtained 
from them may be equally meaningful or even superior as compared to the 
four- dimensional BS theory. 

The authors are indebted to R. A. Brandenburg for his momentum-space 
Faddeev code. This work was supported in part by the NSF-Idaho EP- 
SCoR program under Grant No. Rll-8902065, the NSF Physics program un- 
der Grant No. PHY-8911040, and by the San Diego Supercomputer Center. 
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Table 1: Triton binding energy (in units of MeV) obtained in non-relativistic 
and relativistic calculations using two versions of the relativistic Bonn two- 
nucleon potential. 





Number 






Relativistic 


Potential 


of 
Channels 


Non- Relativistic 


Relativistic 


and 
Charge-Dependent 




5 


8.37 


8.53 


8.35 


Bonn A 












34 


8.32 


8.51 


8.33 




5 


8.16 


8.33 


8.15 


Bonn B 












34 


8.14 


8.34 


8.16 
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Figure 1: Graphical representation of the three-body rearrangement colhsion 
^(j^) ~^ k{ij). Double lines represent interacting pairs. The underlying time 
axis is horizontal, pointing left into the future. 
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